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Abstract 

Self-consistent solutions to the spinor, scalar and BI gravitational field equa- 
tions are obtained. The problems of initial singularity and asymptotically 



^ I isotropization process of the initially anisotropic space-time are studied. It 



is also shown that the introduction of the Cosmological constant (A-term) in 
the Lagrangian generates oscillations of the BI model, which is not the case 
in absence of A term. 
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5^ ; 1. INTRODUCTION 

The problem of initial singularity still remains at the center of modern day cosmology. 
Though Big Bang theory is deep rooted among the scientists dealing with early day cosmol- 
ogy, it is natural to look back if one can model a Universe free from initial singularity. 

The nonlinear generalization of classical field theory remains one of the possible ways 
to over come the difficulties of the theory which considers elementary particles as mathe- 
matical points. The gravitational field equation is nonlinear by nature and the field itself is 
universal and unscreenable. These properties lead to definite physical interest for the proper 
gravitational field to be considered. 

FRW models are widely considered as good approximation of the present and early stages 
of the universe. But the large scale matter distribution in the observable universe, largely 
manifested in the form of discrete structures, does not exhibit homogeneity of a higher 
order. Recent space investigations detect anisotropy in the cosmic microwave background. 
The theoretical arguments Q] and recent experimental datas that support the existence of an 
anisotropic phase that approaches an isotropic one leads to consider the models of universe 
with anisotropic back-ground. Zel'dovich was first to assume that the early isotropization 
of cosmological expanding process can take place as a result of quantum effect of particle 
creation near singularity 0. This assumption was further justified by several authors 
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A Bianchi type-I (BI) Universe, being the straightforward generahzation of the flat 
Freidmann-Robertson- Walker (FRW) Universe, is one of the simplest models of an 
anisotropic Universe that describes a homogenous and spatially flat Universe. Unlike the 
FRW Universe which has the same scale factor for each of the three spatial directions, a 
BI Universe has a different scale factor in each direction, thereby introducing an anisotropy 
to the system. It moreover has the agreeable property that near the singularity it behaves 
like a Kasner Universe, even in the presence of matter, and consequently falls within the 
general analysis of the singularity given by Belinskii et al 0. Also in a Universe filled with 
matter for p = C^? C < 1? it has been shown that any initial anisotropy in a BI Universe 
quickly dies away and a BI Universe eventually evolves into a FRW Universe [0. Since the 
present-day Universe is surprisingly isotropic, this feature of the BI Universe makes it a 
prime candidate for studying the possible effects of an anisotropy in the early Universe on 
present-day observations. 

Nonlinear spinor field (NLSF) in external FRW cosmological gravitational field was first 
studied by Shikin [§ . Here we would like to note that the presence of a singular point to time 
in its space-time metric is another important property of the isotropic model. The presence 
such a singular point means that the time is restricted. The motivation for introduction 
a nonlinear term in the spinor field Lagrangian was to answer the natural question that 
arises in connection with the presence of singular point, i.e., to what extent the presence of 
a singular point is an inherent property of the relativistic cosmological models or whether 
it is only a consequence of specific simplifying assumptions underlying these models. The 
study by Shikin shows that the presence of spinor filed nonlinearity is not enough to 
remove singularity in a FRW space-time. The natural choice was to introduce anisotropy 
in the model and analyze the nonlinear spinor field equations in an external BI universe 
that was carried out in |P,|lO|,|lT],|l2| Recently, we study [|13|,|14| the role of the cosmological 
constant (A) in the Lagrangian which together with Newton's gravitational constant (G) is 
considered as the fundamental constants in Einstein's theory of gravity. We also considered 
an interacting system of spinor and scalar field in a BI space-time |]T3,rB,r7|. 

In this report we consider a self-consistent system of spinor, scalar and Bianchi type-I 
gravitation fields in presence of perfect fluid and A-term in an attepmt to 
Obtain exact solutions to the nonlinear spinor, scalar and gravitational fields; 
Study initial singularity and possibilities of its elimination; 
Isotropization of initially anisotropic cosmological model; 
Possibility of the existence of oscillating models; 
Role of A-term and perfect fluid in the evolution of BI universe. 



2. FUNDAMENTAL EQUATIONS AND GENERAL SOLUTIONS 

The action of the nonlinear spinor, scalar and gravitational fields can be written as 

5(^;V',^,¥^)= / Ly^dn (2.1) 

with 



L — + Lgp + Lsc + Lint + -^m- 



(2.2) 



Here Lgp is the spinor field Lagrangian given by 



— mijjijj. 



(2.3) 



The mass-less scalar field Lagrangian is chosen as 
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Lsc = l^,a^''' (2.4) 
while the interaction Lagrangian we choose |T0|JT5|JT7[] 

i^int = ^¥^,a¥^'"i^, (2.5) 

with A being the coupling constant and F is some arbitrary functions of invariants generated 
from the real bilinear forms of a spinor field. We choose the nonlinear F to be the function of 
I = = {ipipy and J = = [iip'j^tljy , i.e., F = F{I, J), that describes the nonlinearity 
in the most general of its form fl^ . 

Lg corresponds to the gravitational field 

where R is the scalar curvature, n = 8ttG with G being the Einstein's gravitational constant 
and A is the cosmo logical constant. The gravitational field in our case is given by a Bianchi 
type I (BI) metric in the form 

ds'^ = dt^ - a^{t)dx^ - b\t)dy^ - c\t)dz\ (2.7) 

with a{t), b{t), c{t) being the functions of time only. 

Variation of ( |2.1D with respect to spinor field ip {ip) gives nonlinear spinor field equations 

iYV^^ -mil) + Vi) + Qi'^^i) = 0, (2.8a) 

iV^ipY + m^-V^- gi^-f^ = 0, (2.8b) 

where we denote 

OF OF 

whereas, variation of ( |2.1| ) with respect to scalar field yields the following scalar field equation 

^ d_ 

y/—g dx^ 

Varying (|2.1|) with respect to metric tensor g^y one finds the Einstein's field equation 



gg^^\l + XF)^A=Q. (2.9) 
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(2.10a) 


KTi - A, 


(2.10b) 


kTI - A, 


(2.10c) 




(2.10d) 



where point means differentiation with respect to t and Tj^ is the energy-momentum tensor 
of the material field given by 
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rpu rpU:rpU:rp^:rpL, fOll) 

-^11 sp 11 ^ sc fj. ^ mt fj. ^ m fi- K'^-^^J 

Here T^p^ is the energy-momentum tensor of the spinor field 

^sp'^ = la'" (v^TmV.V' + V^7.V^^ - V^^7.^ - - S',L,p (2.12) 

where L^p with respect to (|2.8|) takes the form 

L,p = -{VS + gP). (2.13) 
The energy-momentum tensor of the scalar field T^^fi is given by 

Tse; = - 5;:^sc. (2.14) 

For the interaction field we find 

T;:,^ = XFcp^^cp'^-6;U^,. (2.15) 

'^fl(m) — ~Py ~Pi ~P) is the energy-momentum tensor of a perfect fiuid. Energy e is 
related to the pressure p by the equation of state p = (e. Here ( varies between the 
interval < ( < 1, whereas C = describes the dust Universe, C = | presents radiation 
Universe, | < < 1 ascribes hard Universe and ( = 1 corresponds to the stiff matter. 
The explicit form of the covariant derivative of spinor is IT^JH 



V,^ = ^-r,^, (2.16a) 

V,^ = ^ + V^r„ (2.16b) 

where r^(a;) are spinor affine connection matrices. For the metric ( p.7|) one has the following 
components of the affine spinor connections 

ro = o, ri = id(t)7V, T, = h{t)ff, Ts = ^cit)ff. (2.17) 



We study the space-independent solutions to the spinor and scalar field equations ( pl8| 
and (|2.9|) so that ip = ip{t) and ip = (p(t). defining 



r = a6c=v^ (2.18) 
from (|2.9|) for the scalar field we have 

The spinor field equation ( p. 8a ) we rewrite as 

i'y°(— + —]'^-mip + Vip + gi'y^ip = 0. (2.20) 
Vat 2tJ 

Setting Vj{t) = y/Tipj{t), j = 1,2,3,4, from (|2.2CI| ) one deduces the following system of 
equations: 
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Vi+i{m-V)Vi-gVs = 0, 
V2 + %{m - V)V2 - = 0, 
- i{m - V)Vi + gvi = 0, 
Vi - i{m - V)Vi + gV2 = 0. 



(2.21a) 
(2.21b) 
(2.21c) 
(2.21d) 



Prom ( 2.8a ) we also write the equations for the invariants S, P and A = ipj^j^ilj 

So-2gAo = 0, 
Po - 2(m -V)Ao = 0, 
Ao + 2{m-V)Po + 2gSo = 0, 

where 5*0 = tS, Pq = tP, and Aq = tA, leading to the following relation 

S^ + P^ + A'' = CyT\ 



(2.22a) 
(2.22b) 
(2.22c) 



const. 



(2.23) 



Giving the concrete form of F from (|2.21|) one writes the spinor functions in explicit 
form and using the solutions obtained one can write the components of spinor current: 



The component 



T 



■[v:v, + v*v2 + v*vs + v:v,] 



(2.24) 



(2.25) 



defines the charge density of spinor field that has the following chronometric-invariant form 

^= (Jo ■/)'/'. (2.26) 
The total charge of spinor field is defined as 



Q = J gdxdydz = grV, 



(2.27) 



where V is the volume. From the spin tensor 



one finds chronometric invariant spin tensor 



^ch ~ Wij,0^ 



,0\l/2 



and the projection of the spin vector on k axis 

Sk= J S^'^gdxdydz = Sii;\V. 



:2.28) 



(2.29) 



(2.30) 



Let us now solve the Einstein equations. To do it we first write the expressions for the 
components of the energy-momentum tensor explicitly: 
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mS + 



2r2(l + \F) 



(2.31) 



Tl 



rp2 

^2 



rp3 



vs + gp 



p. 



2r2(l + XF) 

In account of ( p.31|) from (|2.10a|) , (|2.10b| ),( p.l0c|) we find the metric functions |12 



a{t) 

m 

c{t) 



2X1 + X3 



dt 



r{t) 



{DiDl 



Xi — X3 



Xi + 2X3 



dt 

W) 
dt 



r{t) 



(2.32) 
(2.33) 
(2.34) 



As one sees from ( p.32|) , (|2.33| ) and ( p.34|) , for r = t" with n > 1 the exponent tends to unity 
at large t, and the anisotropic model becomes isotropic one. Let us also write the invariants 
of gravitational field. They are the Ricci scalar Ii = R ^ 1/t, I2 = Rf^uR^'^ = R'^R'i ~ 
and the Kretschmann scalar = R^p^yR'^^^'^ ~ As we see, the space-time becomes 

singular at a point where r = 0, as well as the scalar and spinor fields. Thus we see, all the 
functions in question are expressed via r. In what follows, we write the equation for r and 

study it in details. 

Summation of Einstein equations (p.lOa|) , (|2.10b|) , (p.l0c|) and (|2.10d|) multiplied by 3 
gives 



- = ^K{niS + vs + gp 

T 2 



P 



3A. 



:2.35) 



For the right-hand-side of ( p.35| ) to be a function of r only, the solution to this equation is 



well-known 20 



From energy-momentum conservation law 



rpu 



rpu I -pi/ rpp -pp rpu 



in account oi p = C,e we obtain 



r 



p 



(2.36) 



^2.37) 



In our consideration of F as a function of /, J or / ± J we ge t these arguments, as well as 
V and g as functions of r. Hence the right-hand-side of ( p.35|) is a function of r only. Then 
( p.35| ), multiplied by 2f can be written as 

(2.38) 



2f f = [3(K(mS + VS + gP + e-p)- 2A)r]f 
Solution to the equation ( |2.38| ) we write in quadrature 

dr 



r r 



/M/(r)dr 



(2.39) 



From here on we set t p = , as it only gives the shift of the initial time. Given the explicit 
form of F{I , J), from ( |2.39 ) one finds concrete function T{t). Once the value of r is obtained, 
one can get expressions for components ipj{t), j = 1,2,3,4. Thus the initial systems of 
spinor, scalar and gravitational field equations are completely integrated. 
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3. ANALYSIS OF THE RESULTS 

In this section we analyze the system of spinor, scalar and gravitational field for some 
concrete form of F . 

As it was mentioned earlier, we choose F to be a function of / and J only. We first 

consider the case with F = F{I), that will followed by F = F{J). 

case (i) Let us consider the case with F = F{I). In this case from (|2.22a| ) one finds 

S = Cq/t, Co = const. (3.1) 



In this case for spinor field components we find |12 



(3.2) 



with Ci being the integration constants and are related to Co as Cq = Cl + C"^ — — Cl- 

Here (3 = /(m - V)dt. 

For the components of the spin current from ( 2.24 ) we find 

f = l[Cf + C| + Ci + Cj], = —[CiC + C2C3]cos(2/5), 
r ar 

f = |-[CiC4 - C2C3]sin(2/?), f = -[CiCs - C2C4]cos(2/3), 
or CT 

whereas, for the projection of spin vectors on the X, Y and Z axis we find 

^23,0 ^ '-"l'-"2 + '-^3'-^4 ^31,0 ^ g ^12,0 ^ '-'1 ~ '-"2 + '-"3 ~ '-"4 

hcT ' ' 2a6r 

Total charge of the system in a volume V in this case is 

Q = [Cl + Cl + Cl + Cl]V. (3.3) 

Thus, for r 7^ the components of spin current and the projection of spin vectors are 
singularity-free and the total charge of the system in a finite volume is always finite. 
Let us set F = S"'. In this case for r we write the second order equation 

f = ^(m + Anr"-V2(A + r")2 + £o(l-C)M) -3Ar, (3.4) 



or the one of first order 

dr 



K[mT + Ar"/2(A + r")] - Ar^ + eot^''^ + yf 



V3t. (3.5) 



Here we set C = 1 and Co = 1. It should be emphasized that, r cannot be negative, i.e., 
it has a lower limit. On the other hand, for A > 0, for the integrant of ( p.5| ) to be positive 
definite, r should have some upper limit, yi in ( |3.5| ) is an arbitrary integration constant. 
Giving boundary conditions the equation (|3^) can be solved by continuous analog of Newton 
method |2l|, whereas given initial (or asymptotical) value one can numerically solve ( p.5|) 
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by Runge-Kutta method. Since yi in ( p.5|) is arbitrary, we use continuous analog of Newton 
method. 

Let us formulate the boundary condition. At initial stage perfect fluid plays the principal 
role and r at this stage obeys 



3k 



^o(i-c)M, 



with the solution 



ri + C 



2/(lK) 



Ci = const. 



whereas, at large t we have for r 



with the solutions 



3ft; 



r 



m — 3Ar, 



(3.6) 



(3.7) 



(3.8) 



B 



(3/4)ft:mt2 



A = 0, 



sin(V3At) + (ft;m/2A), A > 0, 



(3.9) 



sinh(V^3At) + (ft:m/2A), A < 0. 



Note that, for the solutions to be regular, r should be non-zero. On the other hand, being 
\f--gi T is non-negative. As one sees, in this case the right hand side of ( |3.4| ) is a continuous 
function. A positive c\ in ( |3.7] ) provides us non-zero r as t — > 0, whereas, manipulating the 
parameters ft;, m, A in ( |3.9| ) we always have a positive r as t — oo. 
Further, setting 



r = e + [(5 - A)l{t^ - to)]t + {InA - tQB)/{tN - to) 
into (|3.4| ) we get the second order equation for ^{t) 



(3.10) 



(3.11) 



where Phi{^) coincides with the right hand side of ( p.4| ) with r substituted by ( p.lO| ). The 
boundary conditions in this case regardless to the value of A read 



m\ 



t = tQ- 



0, m\ 



t=tr 



0. 



Let us now solve ( p.ll| ) together with (|3.12|) . For this we rewrite (|3.11| ) in the form 



^(0 



m = 0. 



Assume 

^ = ^{t,a), 0<a<+oo, 
with a being some parameter. Further we demand 



(3.12) 



(3.13) 



(3.14) 
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da 



In this case we obtain 



nat,cr)] = nat,o)]e- 



(3.15) 



(3.16) 



It means for a large enough we tend to our initial pro blem. Let us assume that ^ possesses 
continuous differentiation. Then the left hand side of (|3.15| ) can be written as 



where we define 



da 

Thus we now have the following equation for t]: 

d'^r] a$ ,d^^ ^^^^^ 



with the boundary conditions 



r]{0,a) = 0, r]{oo,a) 



0. 



(3.17) 



(3.18) 



(3.19) 



(3.20) 



Thus we now have the boundary value problem ( |3.19|) and ( |3.20| ). To solve ( p.l9| ) we need 
to give some initial value of ^. Setting 



^{t,a = 0)=Co, 
from (|3.19|) we find ri{t, 0). Further from ( p.l8|) we find 

e(t,cT) =eo + (Aa)r7(t,0), 



(3.21) 



(3.22) 



where Aa is the step along a and lies between [0, 1], i.e., < Aa < 1. For Act = 1 we have 
the classical Newton method. This process is performed again and again until rj becomes 
small enough. 

Function is completely problem depended. Since ^ = at both left and right boundary, 
we have to choose that is zero at the end-points. In our particular initial 
approximation of ^ one may consider any non-negative function that tends to zero at to 
{t = 0) and tjy {t = oo), e.g.. 



e(t) = (t-to) X {tN-t), 



or, 



TT 



sm 



't]\f — to 



t). 



(3.23) 



(3.24) 



Note that the equation ( ^.4] ) with boundary conditions ( |3.7| ) and ( |3.9| ) is a multi-parametric 
problem. The parameters in questions are the Einstein s gravitational constant k, spinor 
mass m, coupling constant A, power of nonlinearity n, (, cosmological constant A and the 
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constant of integration Ci. We choose k = 1 and Ci = l.O-E" — 10 for all cases. It should 
be emphasized that the power of nonlinearity n may be any number including a negative 
one. In our case we set n = 2. Note that r bahaves more or less in the same manner when 
n is even, whereas, for odd n we have totally different picture depending on its concrete 
value. For spinor mass we put m = 1. The coupling constant A may be both positive and 
negative and plays a crucial role here. We set A = 1 and A = for interacting and case 
with minimal coupling, respectively. The parameter ( describes perfect fluid and takes the 
values C = (dust), C = 1/3 (radiation), ( = 2/3 (hard universe) and C = 1 (stiff matter). 
In the figures the notations "d", "r", "h" and "s" stand for dust, radiation, hard universe 
and stiff matter, respectively. Finally, for cosmological constant we set A = (—0.4,0.0,0.4). 
Changing these parameters we can manipulate the evolution of r, e.g., for A > we can 
obtain r that is positive everywhere. The numerical experiment with different values of 
these parameters we leave for some future papers. Here we solve the problem for some 
definite values mentioned earlier and chalk out the main feature. For a positive A we always 
have oscillatory mode, whereas for negative A we have inflation-like solution. In absence of 
A-term we have the solution that steadily approaches to isotropic stage for a massive spinor, 
while as will be shown later in case of F = P", for a massless spinor field isotropization 
does not take place. It should be mentioned that interaction term plays no significant role 
in this case, hence we have the identical picture for minimal coupling and direct interaction 
between spinor and scalar fields. Spinor mass plays a crucial role in absence of A-term. Note 
that, for a negative A-term evolution of r is completely dominated by A, and independent 
to the value of m, n and ( we always have the same picture. As it was mentioned earlier, 
at a space-time point where r = 0, we have singularity, whereas, in the intervals where 
r < 0, corresponding solution is not a physical one. We would like to remark that, giving 
suitable values for the parameters, it is possible to find solutions those are everywhere non- 
negative, but singularity-free solutions can be achieved only for A > 0. It should also be 
emphasized that, regular solutions may be obtained by virtue of spinor field nonlinearity |12 



or interaction term |jT^ , but this singularity-free solutions are attained at the cost of broken 



dominant energy condition in Hawking-Penrose theorem p^j23|l . 



The energy density of the system in this case has the form 

,0 mCo , CV"-^ , So 



r +2(r" + AQ) + ri+C- ^^'^^^ 



For the total energy we then have 



E = I T^^-^gdxdydz = [mCo + ^ + V. (3.26) 



As one sees from ( |3.26| ), if r is nontrivial, then the total energy of the system in a finite 
volume V is always finite. 

case (ii) Let us consider the case when F = F{J). In the case considered we assume 
the spinor field to be massless. It gives T> = 0. Let us note that, in the unified nonlinear 
spinor theory of Heisenberg, the massive term remains absent, and according to Heisenberg, 



the particle mass should be obtained as a result of quantization of spinor prematter ||2^. In 
the nonlinear generalization of classical field equations, the massive term does not possess 
the significance that it possesses in the linear one, as it by no means defines total energy 
(or mass) of the nonlinear field system. Thus without losing the generality we can consider 
mass-less spinor field putting m = 0. Then from ( |2.22bD one gets 

P = Do/t, Do = const. (3.27) 



In this case the spinor field components take the form 
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"03 



V'4 



The integration constants Di are connected to Dq by Do = 2 {D\ + Z}| — Z}| — D 

set cr = / Qdt. 

For the components of the spin current from (|2.24|) we find 

f = -[Dl + Dl + Dl + D% j' = —[D^D^ + DiD4]cos(2a), 
r ar 

j2 = -i[D2A3 - Z^iD4]sin(2or), f = —[DiD; - D2D^]cos{2a), 

OT CT 

whereas, for the projection of spin vectors on the X, Y and Z axis we find 

2(^1^2 + D3D4) ,31,0 n 012,0 Dl-D-^ 



S 



23,0 



S' 



0, s 



J^4 



(3.28) 



Here we 



bcT 2abT 
Let us now choose F = P". In this case in view of massless spinor field for r we obtain 



f = ^(Anr"-V2(A + r")^ + 5o(l - O/r^ 



3Ar. 



At t = we have the boundary condition as in the previous case, i.e.. 



1 + C 



2/(1+0 



Ci = const. 



(3.29) 



(3.30) 



whereas, at large t we have 



it) 



qit + q2, A = 0, 
sin(V3At), A > 



sinh(V^^t), A < 0. 



(3.31) 



Here gi, q2 are constants of integration. We solve the equations (|3.29|) , (|3.30| ) and (p.31|) 
numerically using continuous analog of Newton method. 

It should be mentioned that the results for C = and C = 1 obtained with |?7| < 
l.OE — 13. Iteration in this case varies from 20 to 80. For ( = 1/3 and C = 2/3 the 
value of |?7| is rather rough. In case of massless spinor field within the scope of parameters 
chosen here, no solutions were available with a positive A for radiation and hard universe. 
Corresponding solutions are presented graphically in the figures followed. As one observes 
in the following figures, solutions are singular at initial time and isotropization process takes 
place for massive spinor field. For positive A-term solutions are always oscillatory and the 
amplitude of oscillations decreases as the value of ( increases. In absence of A-term beside 
the spinor mass perfect fluid plays a significant role. Especially, for dust and radiation, 
expansion of r starts beginning from some positive value of t. 
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FIGURES 




FIG. 1. Evolution of r in case of minimal coupling with positive A. 




FIG. 2. Evolution of r for a positive A term with direct interaction between spinor and scalar 
fields. As one sees, interacting term is significant only at the vicinity of t = 0. 




FIG. 3. Evolution of r for a positive A term in absence of spinor mass. Only two solution are 
available (dust and stiff matter). 
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FIG. 4. Evolution of r absence of A term. For dust and radiation expansion starts at some 
positive value of t, and the rate of expansion is rather rapid. 



m 




FIG. 5. Evolution of r in absence of A and massive terms. There is no isotropization in the 
case concerned. 




FIG. 6. Evolution of r with a negative A. Effects of mass, interaction and perfect fluid are 
completely suppressed by A-term. 
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FIG. 7. Evolution of r with a positive A for different values of n (power of nonlinearity) . As 
one sees, for n even, the solutions more or less coincide, while it is not the case when n is odd. 
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FIG. 8. Evolution of r when BI universe is filled with perfect fluid only. 



4. CONCLUSION 

Within the framework of the simplest model of interacting spinor and scalar fields it is 
shown that the A term plays very important role in BI cosmology. In particular, it invokes 
oscillations in the model which is not the case when A term remains absent. It should be 
noted that regularity of the solutions obtained by virtue of A term, specially for the linear 
spinor field does not violate dominant energy condition, while this is not the case when 
regular solutions are attained by means of nonlinear term. Note that in presence of A-term 
the role of other parameters such as order of nonlinearity n, perfect fluid parameter ^ and 
spinor mass in the evolution process are rather local, while the global process are totally 
determined by the A-term, e.g., for a positive A we have always oscillatory mode, while for 
a negative A solution is always inflation-like. 



Spinor field with induced nonlinearity in a Bianchi type-I Universe 
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